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ABSTRACT 

The color magnitude diagram (CMD) of NGC 1851 presents two subgiant branches 
(SGB), probably due the presence of two populations differing in total CNO content. 
We test the idea that a difference in total CNO may simulate an age difference when 
comparing the CMD of clusters to derive relative ages. We compare NGC 1851 with 
NGC 6121 (M4), a cluster of very similar [Fe/H]. We find that, with a suitable shift of 
the CMDs that brings the two red horizontal branches at the same magnitude level, 
the unevolved main sequence and red giant branch match, but the SGB of NGC 6121 
and its red giant branch "bump" are fainter than in NGC 1851. In particular, the 
SGB of NGC 6121 is even slightly fainter than the the faint SGB in NGC 1851. Both 
these features can be explained if the total CNO in NGC 6121 is larger than that in 
NGC 1851, even if the two clusters are coeval. We conclude by warning that different 
initial C+N+O abundances between two clusters, otherwise similar in metallicity and 
age, may lead to differences in the turnoff morphology that can be easily attributed 
to an age difference. 
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1 INTRODUCTION 

The ages of Globular Clusters (GCs), the oldest systems 
populating our Galaxy for which an age estimate is feasible, 
provide a lower limit to the age of the Universe, and therefore 
constitute one of the important subjects of stellar astrophys- 
ical research. Unfortunately, absolute age measurements are 
still affected by a large number of uncertainties, in particu- 
lar due to significan t uncertainties in GC distances and red- 
dening values (e.g. iGratton et al.l [20031 ). Attention is then 
mostly devoted to the determination of relative ages. These 
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can be obtained with much higher precision, by measuring 
the position of the main-sequence (MS) turnoff (TO) -which 
is the most-used age indicator- relative to some other feature 
in the color-magnitude diagram (CMD) , whose location has 
little or no dependence o n age, like the horizon tal branch 
(HB) (|Stetson et al.lll99rj ; ISaraiedini et ai]|l997l ). Relative 
ages provide fundamental information on the formation and 
early evolution of the Galaxy. 

Several procedures, and better and better data sets have 
been employed to derive relative ages (|Buonanno et al. 1998; 
Salaris fc Weiss! ll99Sl; iRosenberg et all 1 19991; IVandenBergl 
200d : iDe Angeli et al.l 120051; iMarin-Franch et alj l2009af r 
The latest work bv lMarin-Franch et al. (2009a) employs the 
homogeneous photometric database achieved for 64 GCs in 
the ACS Survey of Galacti c GCs (the HST Treasury pro- 
gram [Sarajedinret^LjllQOj)), and compares the relative po- 
sitions of the clusters' MS TOs, using MS fitting to cross- 
compare clusters within the sample. According to the au- 
thors, the method provides relative ages to a formal preci- 
sion of 2% - 7%. One of the problems of the Marin- Franch 
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et al. study is that it does not include any consideration of 
possible systematic elemental abundance differences among 
GCs that apparently have the same [Fe/H]. This problem is 
clearly to be attacked, now that we know that GCs may 
harbour different stellar populati ons, whether because of 
differences in the helium content llD'Antona fc Caloil kotM; 
lNorrisll20o4 iD'Antona et al.ll2005l ; IPiotto et al.ll2005l.l2007f); 
differ ences in a-element abundances in general (|lvans et al.l 
1999); or differences specifically in total C+N+O, as will 
be discussed below. Limited helium content differences (in 
mass fraction, SY ~ 0.04 — 0.05) do not seem to affec t 
the TO luminosity at a given age l|D'Antona et alj|2002h . 
and so in principle do not affect the relative age determina- 
tion. In contrast, stellar modelling shows that an increase 
in the total C+N+ O content shifts the isochrone turnoffs 
to lower luminosity ( Bazzano et al ] |l982l ;l Cassisi et al . 2008; 
I Ventura et al.ll2~009a ) . without significantly changing the po- 
sition of the unevolved MS, giant branch or HB. 

Recently we have obtained observational evidence of 
the influence of such an additional parameter. In fact, 
two well separated subgiant branches (SGB) are presen t 
in the CMD of the cluster NGC 1851 (|Milone et alj|2008h . 
The faint SGB may be interpreted as being populated by 
stars either older by ~lGyr, or having equal (or slightly 
youn ger) age and an increase in the total C +N+O con- 
tent ijcassisi et all 120081 ; IVentura et~aH l2009al ). Formation 
models for multip le populations in GCs are still very much 
in their infancy dBekki fc Norrisl 120061 ; iBekki et all 120071 ; 
iD'Ercole et al.ll2008l ). but no current model can easily ac- 
comodate a second stellar generation younger by 1 Gyr. On 
the other hand, a total C+N+O spread among the stars 
in NGC 1851 does have a n independent ob servational basis 
l|Yong fc Grundah]||2008l ; lYong et al.l 120091 ') . An interpreta- 
tion based on a different total C+N+O content is a variant of 
the "two stellar generations" models for GC formation that 
have already been investigated for the case of multiple main 
sequences, possibly showing differences in helium content. In 
this scheme, the second stellar generation forms out of mat- 
ter polluted by the stellar winds of m assive asymptotic gi- 
ant branch stars of the first generation |Cottrell fc Da Costal 
I 19811 ; IVentura et al.ll200ll ). althou gh many details of su ch a 
model still need to be worked out (|D'Ercole et al.ll2008T ). 

Now that we have such observational evidence, we can 
pose the question whether relative age indicators already 
present in the literature can be interpreted differently. In 
the following, we explore the case of the relative ages of 
the two clusters NGC 1851 and NGC 6121 (M4), that have 
always been considere d twins from the point of view of 
meta l abundance (e.g. Zinn fc Westlll984l; iRutledge et all 
Il997l ; iRosenberg et al.1 Il999l ; iMarin-Franch et al.l 1200981 ). 
with NGC 1851 somewhat younger if age is the only sec- 
ond parameter. 

The influence of varying the CNO abundances in stel- 
lar models is well documented in the astrophysical litera- 
ture (e ; g_ ; _SJjn£d^_&_ftenJ_970, Bazzano et al. 1982; see 
also IVandenBerg fc Belli (|200ll ) and references therein) . In 
the last decade, however, the attention was mainly devoted 
to the effect of variations in oxygen -and in the other a- 
elements, following the discovery th at population II stars 
had a larger than solar [O/Fe] (e.g. lMcWilliamlll997l) . An 
increase by 0.3 dex in [O/Fe] implies a 1 Gyr reduction in 
age for a fixed turnoff luminosity. Recently the attention has 
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Figure 1. Full view of the CMD diagrams in V vs. B-I of the 
clusters NGC 1851 (top figure) and M4 (bottom), superposed 
to the isochrones by applying magnitude and color shifts as ex- 
plained in the text. For age t=12 Gyr and Y=0.24, the standard 
(CNOxl) Z=0.0010 isochrone (green) and the isochrone CNOx5 
(Z=0.0035, yellow) are shown as full lines. The mixtures are de- 
scribed in Ventura et al. (2009a). The full (blue) line is the ZAHB 
locus of models by Ventura et al. (2009b) for Z=10 -3 and Y=0.24, 
the dash— dotted (red) line is the (slightly dimmer) ZAHB locus 
for the Y=0.24, CNOx5 mixture. The black lines on the blue 
side of the HB are the ZAHB locations for the mixture CNOx3, 
Y=0.26 (dashed) and CNOx5, Y=0.28 (dash-dotted). 



been focused on the possible differences in total C+N+O de- 
riving from the presence of multiple stellar generations, so 
different mix t ures of these elements have been considered. 
ICassisi et al.l l|2008T ) adopt an enhanced mixture based on a 
combination der ived from some obser vations of abundance 
anomalies, while lVentura et al.l (|2009al ) resorted more on the 
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Table 1. Relative ages for NGC 1851 and NGC 6121 



Name 


[Fc/H] CG 


relative ages 


relative ages 






Dc Angeli 2005 


Marin Franch 2009 


NGC 1851 


-1.03 


0.82 ± 0.07 


0.78 ± 0.04 


NGC 6121 


-1.05 


0.91 ± 0.04 


0.98 ± 0.05 
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chemical yields provided by theoretical massive AGB mod- 
els. In any case, we are still far from a complete investigation 
of the relative role of the three elements, that may become 
more necessary following the results of the present work. 
We base our considerations mainly on the models computed 
to explain the double subgiant branch in NGC 1851. This 
problem is ver y timely, as it is also confir med by the very 
recent exam bv lMarin-Franch et al.l (|2009bl ) of the influence, 
of CNO and helium variations, on the relative age determi- 
nation of GCs, purely from a theoretical point of view. 
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2 THE TWO CLUSTERS NGC 1851 AND 
NGC 6121 

In this analysis we exploit the archive of hom ogeneous CCD 
photometry maintained by one of us (see, e.g. JStetson et all 
1 19981 ; IStetson~ll200Cl . 12009 ). Th e photometry o f NG C 1851 
has already been described in iMilone et al. I (I2009T I; their 
CMD is reproduced here as the blue points in Fig. 1. Note 
that the ground-based photometry confirms the split of the 
SGB into faint and bright branch es, as originally found 
from HST data (IMilone et al.1 120081 '). In addition, the fit by 
IVentura et al.l |2009al ) of the bright SGB with a standard 
isochrone of 12 Gyr, and the fit of the corresponding faint 
SGB with the same isochrone and C+N+O abundances in- 
creased by a factor three (CNOx3 mixture) is consistent 
with the observations in the V vs. B-I plane (see Fig[2]). 

The photometry of NGC 6121 to be discussed here is 
similarly obtained from analysis of 1,199 individual CCD 
images resulting from 814 exposures (55 of the exposures 
had been taken with the eight-CCD WFI camera on the 
ESO/MPE 2.2m telescope) made during ten observing sea- 
sons on six telescopes (ESO/Dutch 0.9m, CTIO 0.9m, Ja- 
cobus Kapteyn lm, ESO/Danish 1.5m, and Nordic Op- 
tical 2.6m Telescope, in addition to the aforementioned 
ESO/MPE 2.2m). Any given star may have been measured 
in as many as 255 B-, 202 V-, and 253 /-band images. 
The results have been transformed to the photometric sys- 
tem of Lando l t (1992) followi n g me t hodol ogy described in 
IStetson et ail (|l998D ; IStetson1 (|2000l . l2005l ). 

Determination of relative ages requires complex proce- 
dures for determining observational parameters and for com- 
parison with theoretical isochrones. Let us look at the final 
compariso n between the two cluste rs un der discussion in the 
analysi s of |Pe Angeli et all (|2005l ) and iMarm- Franch et al.l 
l|2009al ) summarized in Table 1. Both works make a very 
careful analysis, whose final results do not depend sub- 
stantially on the adopted abun dance scale and isochrones . 
Here we report the result for the lCarretta fc Grattonl (Il997l ) 
abundances, and f or the model s by Cassisi et al.l ( 20041 ) in 
De Angeli, and bv lDotter et~aH (|2007F " in Marfn-Franch. 

The table shows that — according to the age in- 



Figure 2. Location of the TO and MS of the two clusters, for 
the magnitude and color shifts adopted for Fig[T] The standard 
(CNOxl) isochrones for t=12 and 14 Gyr with Z=10~ 3 and 
Y=0.24 are shown as full (yellow) lines. The long-dashed (cyan) 
curve represents the 12 Gyr isochrone for models having Z=10 — 3 
and CNOx3, which fits the faint SGB of NGC 1851 well, while 
the dot-dashed (black) curve is the 12 Gyr isochrone for CNOx5. 



terpretation of the observed differences — NGC 1851 is 
~10% younger than NGC 6121 in the |Pe Angeli et all 
(2005 ) analysis, and ^20 % younger in the analysis by 
iMarin-Franch et al.l (|2009al ). Assuming an age of 12 Gyr for 
NGC 6121, NGC 1851 turns out to be a bit less than 11 Gyr 
(De Angeli), or 9.5 Gyr (Marin- Franch) old, an age differ- 
ence that is significant for models of the formation of the 

Galaxy. 

In their Figure 7, |Pe Angeli et al.1 l|2005h compare the 
V vs. V-I CMDs of these two clusters, showing that, when 
the turnoffs are superposed, the HB of NGC 6121 turns 
out more luminous than the HB of NGC 1851, indicat- 
ing indeed that NGC 6121 must be older. We make the 
reverse procedure. We use the metallicity information, to 
superpose the red parts of the HB: even in the presence 
of multiple populations, the red HB clumps are supposed 
to reflect the composition of the first stellar generation, so 
that their location is presumed to be unaffected, e.g., by 
altered helium. The match is made in the V versus B-I 
plane and is shown in Fig. Q] and in the enlarged view of 
Fig(2] We use for this comparison the absolute magnitude 
M v , by positioning the red clumps (RHB) of NGC 1851 
and NGC 6121 on the zer o-age horizontal bran ch (ZAHB) 
of our theoretical models (I Ventura et al.l [200 9b) . drawn in 
the figure. We selected by eye all the stars that, on the ba- 
sis of their position in the CMD, belong very likely to the 
RHB and calculated their median V magnitude. We find 
Vhb=16.19±0.01 and Vhb=13.44±0.01 for NGC 1851 and 
NGC 6121 respectively. The error associated to the HB po- 
sition is calculated as <j/y^{N — 1) where N is the num- 
ber of red HB stars and a is the 68.26 th percentile of the 
absolute deviation from the average of their V distribu- 
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Figure 3. Differential (top panels in each figure) and cumulative 
luminosity functions of the RGBs in NGC 1851 and NGC 6121. 
For each cluster, we plot two histograms. The black ones are cal- 
culated by counting the red giants into intervals of 0.13 mag in 
V. Blue histograms are constructed along the RGB with a step 
of 0.013 magnitudes (a tenth of the bin size) in V, by shifting 
the bin of fixed width of 0.13 mag, with successive steps of 0.013 
magnitude. The (red) arrows mark the location of the bump. For 
NGC 6121 the bump is clearly identified in both the differen- 
tial and cumulative luminosity functions, while in NGC 1851 the 
bump feature is broader in magnitude. This more complex bump 
structure may indeed reflect the fact that NGC 1851 itself con- 
tains two populations with different C+N+O abundance. 



tion. For the aims of this paper however we do not need 
the individual distance moduli of the two clusters (e.g., we 
could have chosen smaller values, assuming a smaller age) 
but their relative displacement, to which we can then as- 
sociate an error of ±0.01 mag. We also choose color shifts 
so that the red giant branch (RGB) locations coincide. The 



apparent V magnitudes of NGC 1851 are shifted by 6V=— 
15.52mag and the colors are shifted by S(B-I)— -0.08. The 
corresponding shifts for NGC 6121 are 51/=-12.77mag and 
5(B-I)=— 1. On the blue side of the HB, the stars in M4 are 
slightly more luminous than in NGC 1851. This feature is 
not important for the present discussion, but remember that 
a slightly larger helium abundance in the blue HB of M4 may 
lead to greater blue HB luminosity, as shown by the (blue) 
ZAHBs with Y=0.26 (dashed) and Y=0.28 (dash-dotted) 
shown in Fig[T] We see in the figure that the t wo MSs are well 
matche d. The isochrones of 12 and 14 Gyr by I Ventura et all 
l|2009al ) for Z=10~ 3 and standard C+N+O abundances are 
also shown. The TO locations indicate an age diffe r ence o f 
~2 Gyr, as also concluded by iMarm- Franch et all (|2009af) 
Since the SGB and TO of NGC 6121 appear fainter than 
the faint SGB of NGC 1851, another interpretation is possi- 
ble: if NGC 6121 stars have total C+N+O abundance even 
larger than that of the stars on the faint SGB of NGC 1851, 
they could be coeval with the stars in NGC 1 851. This is il- 
lustrate d in Fig [2] by the 12 Gyr isochrone by I Ventura et all 
(2009a) for a composition in which the total CNO is en- 
hanced by a factor of five (CNOx5). 

Are there any independent indications that this solu- 
tion (same age and large difference in total C+N+O abun- 
dances) is to be preferred? A possible hint comes from 
Fig. |31 which shows the differential and cumulative red gi- 
ant luminosity functions, and the location of the red gi- 
ant "bump" due to the penetration of the hydrogen burn- 
ing shell into the mean molecular weight discontinuity (or 
/i-barrier) left by the maximum extension in mass frac- 
tion reached by the convective envelope at the beginning 
of giant-branch evolution. Standard theory and observa- 
tions do not fully agree on the locati on in magnitude of the 
bump. Many models — including the IVentura et al.l (|2009a[ ) 
models — provi de bumps ~0.25mag mo re luminous than the 
observed ones (IZo ccali fc Piottdl200Ch . This may be due to 
uncer tainties either in the metallicity scale (|Zoccali et al.l 
Il999l. e.g.) or in the mixing below the convective envelope 
dGirardi et al ]|200Ch. or t o additional mixing caused by rota- 
tion ijPalacios et al.ll2006l ) . Whatever the causes of this small 
discrepancy, here we are interested mainly in the relative po- 
sitions of the bumps in the two clusters under discussion. Ap- 
plying the distance moduli used to build up Fig[T] to which 
we can assign a relative error not larger than ±0.01 mag, 
we see that the absolute magnitudes (median value) of the 
RG bump are M„=0.65± 0.01 mag and M„=0.79±0.01 mag, 
respectively for NGC 1851 and NGC 6121. To determine the 
bump locations, we selected by hand stars that, according 
to their position in the CMD, have a high probability to 
belong to the RGB bump and calculated their median My 
(My bump ). The error is calculated as: a/ yJ~(N— 1), where a is 
the 68. 26 th percentile of the absolute deviation from the av- 
erage of the My distribution and N is the number of selected 
stars. We see that there is a difference of 0.14±0.02 mag in 
the bump level. Can we change the relative distance of the 
clusters in order to match the bumps? The red HB luminos- 
ity is indeed the best standard parameter for stars having 
the same [Fe/H] . Even if we assume that the HB red clump is 
populated by stars with s imilar [Fe/H] c o ntent b ut different 
C+N+O, the models by IVentura et al] l|2009bl ) show that 
the ZAHB luminosity at the red side of the HB decreases 
only by ~ 0.04 mag, for a five-fold increase in C+N+O, so 
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the HB location can not explain this difference in the bump 
magnitude. An age increase of ~2 Gyr can not explain the 
difference either: according to our models, the bump location 
becomes d immer by only ~0.05 m ag. On the other hand, the 
models by IVentura et al.l (1200981 ) for Z=10~ 3 and different 
C+N+O contents show that the bump luminosity decreases 
(almost) linearly with the C+N+O increase, ft was already 
known that an increase in th e oxygen abundance produ ces a 
lower bump luminosity (e.g. VandenB erg fc Belli 1200 J ). We 
obtain 5Af„(bump)~ 0.037 times the C+N+O abundance 
increase with respect to the "standard", a-enhanced, com- 
position (CNO x f). A difference of 0.14+0.02 mag, for a fixed 
age and helium abundance, corresponds to the difference be- 
tween CNOxl and CNO x (3.8 ± 0.6) bumps. We conclude 
that the difference in absolute magnitude of the bumps is 
consistent with the hypothesis that M4 stars have a C+N+O 
content a factor of ~4 larger than the bright SGB stars of 
NGC 1851. 

Note also that the bump peak is more dispersed in 
NGC 1851 than in NGC 6752. This difference may indeed re- 
flect the fact that NGC 1851 itself contains two populations 
with different C+N+O abundances. 



3 DISCUSSION: A PLEA FOR CNO 

ABUNDANCE DETERMINATIONS AND 
FOR A WIDER THEORETICAL EXAM 

The comparison in Figs. [1] and [2] has shown that the two 
clusters we are considering may differ in age by about 2 Gyr 
(as found bv lMarfn-Franch et al. 2009a). However, now that 
we know that NGC 1851 harbours two SGBs, and that the 
faint SGB is very probably made up of coeval stars richer 
in total C+N+O than the bright SGB stars, it is tempting 
to propose that the stars of NGC 6121 are also coeval with 
"all" the stars of NGC 1851, but they have a global C+N+O 
abundance larger even than the stars on the fainter SGB in 
the latter cluster. An analysis of the magnitude difference of 
the RGB bump locations in the two clusters, and the rela- 
tive location of their SGBs are consistent with the hypothe- 
sis that the global C+N+O in M4 stars is about a factor ~4 
larger than in the main population of NGC 1851. At present, 
this conclusion must be taken more as a working hypothe- 
sis for further discussions than as a proven fact. We need 
both careful, extensive and homogeneous abundance deter- 
minations of the C+N+O abundances in most GCs, and 
further theoretical analysis of the role of the three different 
elements — C, N and O — in shaping the clust er isochrones. 
Preliminary results of IVentura et alj l|2009ah suggest that 
the most important role is played by the late phases of hy- 
drogen burning close to the turnoff , where hydrogen burning 
shifts from p-p to the CN chain. As well known, however, 
also the variations in oxygen, although they do not affect at 
all the turnoff evolution, where the ON cycle does not play 
a role, affect the opacities and the evo lution by producing 
smaller ages l| VandenBerg fc Belli lioOll ). 

From an observational point of view, it is not unlikely 
that detailed elemental abundances may differ from clus- 
ter to cluster, even for t he same iron abun d ance and global 
a-element abundances. Ilvans et al. I (| 19991 , l200ll ). in their 
study of a large sample of giants in M4 (NGC 6121) and 
M5 (clusters with very similar metallicity), show that on 



average oxygen, silicon, aluminum, barium, and lanthanum 
are overabundant in M4 with respect to M5. The C+N+O 
abundance sum is constant to within the obser vational er- 
rors. T hese results are confirmed and enlarged bv lYong et alj 
(2008a), who find that silicon, copper, zinc, and all s-process 
elements ar e approxima t ely 0.3 dex overabundant in M4 rel- 
ative to M5. lYong et~ail (|2008b) compared the rubidium and 
lead abundances in M4 and M5, finding that M4 has global 
abundances of both Rb and Pb much larger than M5, and 
that there is no star-to-star variatio n in the abundanc es of 
these elements. As mentioned before. IVong et all (I2009T ) find 
a large spread in C+N+O among 4 giants of NGC 1851, con- 
firming independently the interpretation of the two SGBs 
in terms of different CNO abundances. In their Fig. 4, they 
compare C+N+O abundances i n NGC 1851 with those of 
the giants in M4 determined by ISmith et~aH (|2005l '). These 
latter abundances show a little spread, and fall in the mid- 
dle of the distribution of NGC 1851 giants. Although this 
result goes in the direction we suggest, the location of M4's 
SGB relative to the NGC 1851 SGBs would require the CNO 
abundances in M4 to be similar to or greater than the largest 
abundances of NGC 1851. The small sample examined, and 
the possible syst e matic differe nces in the a b unda nce anal- 
ysis [Yong et alj (|2009T ) versus ISmith et all (120051 )] do not 
allow us to derive more stringent conclusions. 

In summary, the opening of Pandora's box of the pe- 
culiar morphologies of so many CMDs of GCs may reopen 
the problem of relative age determination. We should care- 
fully examine the abundances of C+N+O in the clusters, 
especially those showing different ages at similar metallic- 
ities, before the derived age spreads may be attributed to 
the modalities of formation of the Galaxy. 
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